Precision Measurement of the Weak Mixing Angle in M0ller Scattering 
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We report on a precision measurement of the parity-violating asymmetry in fixed target electron- 
electron (M0ller) scattering: Apv = ( — 131 ± 14 (stat.) ± 10 (syst.)) x 10" 9 , leading to the deter- 
mination of the weak mixing angle sin 2 9^ = 0.2397 ± 0.0010 (stat.) ± 0.0008 (syst.), evaluated 
at Q 2 — 0.026 GeV 2 . Combining this result with the measurements of sin 2 9^ at the Z° pole, 
the running of the weak mixing angle is observed with over 6<r significance. The measurement sets 
constraints on new physics effects at the TeV scale. 

PACS numbers: 11.30.Er, 12.15.Lk, 12.15.Mm, 13.66.Lm, 13.88.+e, 14.60.Cd 



Precision measurements of weak neutral current pro- 
cesses at low energies rigorously test the Standard Model 
of electroweak interactions. Such measurements are sen- 
sitive to new physics effects at TeV energies, and are 
complementary to searches at high energy colliders. 

One class of low-energy electroweak measurements 
involves scattering of longitudinally polarized elec- 
trons from unpolarized targets, allowing for the de- 
termination of a parity- violating asymmetry Apy = 
{<Jr — <Jl) / {&r + where ctr(l) is the cross section 
for incident right (left)-handed electrons. Apy arises from 
the interference of the weak and electromagnetic ampli- 
tudes [lj and is sensitive to the electroweak coupling con- 
stants and thus the weak mixing angle #w- 

The electroweak parameter sin 2 0^ is defined as the 
ratio of the electromagnetic to the weak isospin coupling 
constants . Possible new physics contributions at very 
high energy scales can be expressed in terms of their im- 
pact on the measured value of sin 2 6^f. Measurements 
at low momentum transfers Q 2 <C Af§ can have sensitiv- 
ity comparable to high energy collider searches for new 
physics provided sin 2 6$ is measured to better than 1%. 

At such a precision, the variation of the coupling con- 
stants with momentum transfer, a fundamental property 
of gauge interactions referred to as running Q , must be 



taken into account. While the running of the electro- 
magnetic and strong coupling constants has been clearly 
established, it has not been unambiguously observed for 
sin 2 6^ so far. The variation of sin 2 9^ from Q 2 w to 
Q 2 = M| is due to higher order amplitudes involving vir- 
tual weak vector bosons and fermions in quantum loops, 
referred to as electroweak radiative corrections |J, |5j . 

To date, the most precise low-energy determinations 
of the weak mixing angle come from studies of parity 
violation in atomic transitions and measurements of 
the neutral current to charge current cross section ratios 
in neutrino- nucleon deep inelastic scattering 0. In this 
Letter, we present a measurement of the weak mixing 
angle in electron-electron (M0ller) scattering, a purely 
leptonic reaction with little theoretical uncertainty. We 
have previously reported the first observation of Apy in 
M0ller scattering [8j. Here, we report on a significantly 
improved measurement of Apy resulting in a precision 
determination of sin 2 8yJ at low momentum transfer. 

At a beam energy of ~ 50 GeV available at End Station 
A at SLAC and a center-of-mass scattering angle of 90°, 
Apy in M0ller scattering is predicted to be ~ 320 parts 
per billion (ppb) at tree level Q. Electroweak radiative 
corrections UOI an(1 the experimental acceptance reduce 
the measured asymmetry by more than 50%. 
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The principal components of the experimental appa- 
ratus are the polarized electron beam, beam diagnostics, 
the liquid hydrogen target, the spectrometer/collimator 
system, and detectors. They are described in our previ- 
ous publications H, El Eli we discuss them briefly here. 

The high-intensity polarized electron beam, delivered 
in ~ 270 ns pulses at the rate of 120 Hz, passes through 
a 1.57 m long cylindrical cell filled with liquid hydrogen 
[l2^ . Scattered particles with 4.4 < 6i a h < 7.5 mr over 
the full range of the azimuth are selected by the magnetic 
spectrometer while the primary beam and forward 
angle photons pass unimpeded to the beam dump. 

Sixty meters downstream of the target, scattered 
M0ller electrons in the range 13-24 GeV form an 
azimuthally-symmetric ring, spatially separated from 
electrons scattered from target protons (ep scattering). 
The charged particle flux is intercepted by the pri- 
mary copper/fused silica fiber sandwich calorimeter. The 
asymmetry is measured by extracting the fractional dif- 
ference in the integrated calorimeter response for incident 
right- and left-handed beam pulses. 

The calorimeter provides both radial and azimuthal 
segmentation. Four radial rings are uniformly covered 
in the azimuth by 10, 20, 20, and 10 photomultipliers, 
respectively. The three inner rings, referred to as the 
Inner, Middle and Outer M0ller rings are predominantly 
sensitive to M0ller scattered electrons. The outermost or 
EP ring intercepts the bulk of the ep flux. 

The background within the M0ller rings is estimated 
to be ~ 8%. It is dominated by radiative ep scatter- 
ing, while neutral particles and charged pions contribute 
less than 1%. Quartz detectors placed behind the M0ller 
detector and shielding record the charged pion asymme- 
try. Target density fluctuations and spurious asymme- 
tries are monitored by intercepting charged particles at 
9\ab ~ 1 rnr with gas ionization chambers [13}. 

The data sample consists of 2.9 x 10 8 and 3.7 x 10 8 
pulses at beam energies of 45.0 and 48.3 GeV respectively, 
collected over three data runs in 2002 and 2003. Roughly 
60% of the data were accumulated in the 2003 run which 
featured a novel "superlattice" photocathode [14| with 
rs 90% beam polarization. 

Data were collected at 120 Hz, with ~ 1 Hz of pulses 
blanked to measure baseline signals. Alternate triggers 
fall into two 60 Hz fixed-phase "time slots" . Within these 
time slots, right-left pulse pairs are formed for indepen- 
dent asymmetry analyses. The helicity sequence in each 
pulse pair was chosen pseudo-randomly. In addition to 
the fast helicity flips, the sign of the electron polarization 
was passively reversed in two independent ways. First, 
the state of a half-wave plate in the laser line was toggled 
each day, guarding against helicity-correlated electronics 
crosstalk. Second, spin precession in the 24.5° bend after 
acceleration created opposite helicity orientations at 45 
and 48 GeV beam energies. For each of the 2002 runs 
the beam energy was changed once, while the change 



was made roughly every four days during the 2003 run. 
Roughly equal statistics were accumulated with oppo- 
site signs of the measured asymmetry, suppressing many 
classes of systematic effects. 

We select pulses with beam charge greater than 10 11 
electrons and require that the beam crosses each beam 
position monitor (BPM) within 1 mm of its geometric 
center. We also require that the beam position and 
charge for each pulse be within 6 standard deviations 
from the running mean value. Typically, the beam charge 
per pulse varies between (4—6) x 10 11 electrons with 0.3% 
pulse-to-pulse jitter, and the beam position is within 
100 jum of each BPM center, with jitter on the order 
of 50 jttm. In order to avoid helicity-dependent biases, 
we reject several pulses before and after a pulse which 
fails a cut. Other than the demand that the beam charge 
asymmetry measured by two independent monitors agree 
to within 10~ 3 , no helicity-dependent cuts are made. 

The right-left asymmetry in the integrated detector 
response for each pulse pair is computed by normalizing 
to incident charge and then correcting for beam fluctua- 
tions. To first order, six correlated parameters describe 
the beam trajectory: charge, energy, and horizontal and 
vertical position and angle. Each parameter is measured 
by two independent monitors, such that device resolution 
and systematic effects can be studied. 

Two methods are used to calibrate the detector sensi- 
tivity to each beam parameter and remove beam-induced 
random and systematic effects from the raw asymmetry. 
One method uses a calibration subset (4%) of the pulses, 
where each beam parameter is modulated periodically 
around its average value by an amount large compared 
to nominal beam fluctuations. The other method applies 
an unbinned least squares linear regression to the pulses 
used for physics analysis. They yield statistically consis- 
tent results to within 3 ppb. Final results are obtained 
with the latter, statistically more powerful technique. 

Additional bias to the measured asymmetry may arise 
from asymmetries in higher order moments of beam dis- 
tributions, such as temporal variations of the beam po- 
sition or energy within a 270 ns beam pulse, coupled to 
the intra-pulse variation of position or energy asymme- 
tries. Such higher order biases are small for the Inner 
and Middle M0ller detector rings, but are observed to be 
significant for the Outer ring. 

During the physics runs, great care was taken to mini- 
mize residual time structure of the beam position at the 
target, keeping it typically within 1 mm. In order to 
measure the possible bias due to such effects, six BPMs 
were instrumented with 23 additional readout channels 
before the 2003 run. Thus, in addition to the average 
beam parameters for each pulse, BPM signals for charge, 
energy, positions, and angles are each digitized in four 
independent time slices (three slices for energy). Correc- 
tions due to intra-pulse variation of beam asymmetries 
are computed by linearly regressing M0ller asymmetries 
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against the beam asymmetries in time slices. 

For the 2003 data, the regression analysis limits the 
possible contribution to the detector asymmetry due to 
the intra-pulse variations to 3 ppb. Since time slice data 
were not available for 2002 data, the Outer M0ller ring 
channels are only used in 2003 data sample. For the 2002 
datasets, these channels set a limit on the maximum pos- 
sible bias in the two innermost M0ller rings (containing 
the bulk of the statistical weight) of less than 5 ppb. 

After linear regression, the integrated responses of all 
the selected rings are averaged to form the raw asym- 
metry A raw . Near-perfect azimuthal symmetry reduces 
the sensitivity to beam fluctuations and right-left beam 
asymmetries. The A raw pulse-pair distribution has an 
average RMS of 215 ppm for the 2002 data and 185 ppm 
for the 2003 data. The cumulative beam asymmetry cor- 
rection is —9.7 ± 1.4 ppb. A correction due to an az- 
imuthal modulation of A raw |l5| from a small non-zero 
transverse component of the beam polarization is found 
to be -3.8 ± 1.5 ppb. 

The average electron beam polarization is P& = 0.89 ± 
0.04, measured every few days by a polarimeter using 
M0ller scattering of beam electrons off a magnetized foil. 
The linearity of the calorimeter response is determined 
to be e = 0.99 ± 0.01 from special calibration runs. 

The physics asymmetry A p h ys is formed from A raw by 
correcting for background contributions, detector linear- 
ity and beam polarization: 



A 



phys 



1 A iaw — AAi 
We 1 - li 



Asymmetry corrections AAi and dilutions fi for various 
background sources are listed in Table [I] 

Figure shows A p h ys for all data, divided into 75 se- 
quential samples in Runs I, II (2002), and III (2003). 
Each Aphys measurement has sign reversals depending 
on the beam energy and the state of the half-wave plate. 
Apv is obtained by correcting each result by the appro- 
priate sign. The combined result is 

A PV = -131 ± 14 (stat.) ± 10 (syst.) ppb . 



TABLE I: Corrections AAi and dilutions /, to A r 
sociated systematic uncertainties. 



and as- 



Source 


AA (ppb) 


f 


Beam (first order) 


-10 ± 1 




Beam (higher order) 


0±3 




Transverse polarization 


-4 ±2 




e~ +p — > e~ +p(+7) 


-7± 1 


0.056 ± 0.007 


e~(-y) +p — * e~ + X 


-22 ±4 


0.009 ± 0.001 


7 + e~ — > e~ ± 7 


Oil 


0.005 ± 0.002 


High energy photons 


3±3 


0.004 ± 0.002 


Synchrotron photons 


Oil 


0.002 ± 0.001 


Pions 


1 ± 1 


0.001 ±0.001 




A pv = -131 ±14 ppb 
X 2 /df = 78.5/74 
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FIG. 1: Aphys for each of 75 data samples. Data collected 
with half-wave plate inserted(removed) at a beam energy of 
45(48) GeV are shown as solid(open) circles(triangles) . The 
solid line represents the grand average, with the expected 
modulation of the asymmetry sign for each beam energy 
and half-wave plate state. Only statistical uncertainties are 
shown. 



In the context of the Standard Model, we interpret 
the measurement of Apy in terms of the effective weak 



mixing angle sin 2 0^(Q) 



A 



PV 



-A{Q\y)p^ [l-4sin 2 ^(Q) + A] 



The average values of the kinematic variables are Q 2 = 
0.026 GeV 2 and y = Q 2 /s ~ 0.6, where s is the square 
of the center-of-mass energy. 



GpQ 2 



V2rra(Q) 1 ± y 4 ± (1 - y) 



f-^QED 



is the effective analyzing power, Gf and a(Q) are the 
Fermi and fine structure constants, respectively |l6| . 

the low-energy ratio of the weak neutral and 
charge current couplings, Tqed = 1.01 ±0.01 is a 
QED radiative correction factor that includes kinemati- 
cally weighted hard initial and final state radiation ef- 
fects and y-dependent contributions from 77 and 7Z 
box and vertex diagrams [171 ]. A contains residual O(a) 
electroweak corrections. The effective analyzing power 
A = 3.25 ± 0.05 ppm is determined from a Monte Carlo 
simulation that accounts for energy losses in the target 
and systematic uncertainties in the spectrometer setup. 

A number of definitions of the low-energy weak mixing 
angle exist [5lll8| and differ in the way various corrections 
of order 0{a) are distributed between terms sin 2 9^J(Q), 
, and A. Here we adopt a definition of the coupling 
w(*9) UM which reproduces the effective leptonic 



epto 

coupling sin 2 &§(M Z ) = sj = 0.23149 ± 0.00015 [3 at 



Z° pole. This implies / e;e ) = 1.0012 ± 0.0005 and A = 
-0.0007 ± 0.0009. We determine at Q 2 = 0.026 GeV 2 

in 2 <7w(Q) = 0.2397 ± 0.0010 (stat.) ± 0.0008 (syst.) 
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FIG. 2: Predicted variation [Tg| of sin 2 6$ as a function of 
momentum transfer Q (solid line) and its estimated theoret- 
ical uncertainty (shaded area). Results of prior low energy 
experiments (g,[l6] (closed triangle, shown at an arbitrarily 
higher Q) and B (op en triangle) are overlaid together with 
the Z° pole value [lal (square) and this measurement (circle). 

Our value is consistent with the Standard Model expec- 
tation 0, [3 sin 2 e$J{Q) = 0.2381 ± 0.0006 and is 6.2a 
away from sin 

&w( M z) (Fig- HI- Interpreting our result 
as a measurement of the electroweak coupling parameter 
sin 2 9 w (M z )ms yields 

sin 2 6' w (M z ) M g = 0.2330 ± 0.0011 (stat.) 

± 0.0009 (syst.) ± 0.0006 (theory). 

The last uncertainty is from the evolution to Mz- 

Our measurement of Apy can also be used to set 
limits on the size of possible new contributions beyond 
the Standard Model. Quite generally, we set a limit 
on the scale All of a new left-handed contact inter- 
action characterized by a term in the Lagrangian [lflj l 
C = ±(47r/2A L r2 )(eL7 /i e L ). At 95% C.L. a tree-level cal- 
culation yields A+ L > 7 TeV and A^l > 16 TeV, for 
potential positive and negative deviations, respectively. 
As an example of a specific model, the 95% C.L. on the 
mass of Z x boson appearing in the grand unified model 
with SO(10) symmetry HEll is M Zx > 1.0 TeV. 

In summary, we have reported a new measurement of 
Apy in M0ller scattering with an accuracy of 17 ppb. 
This leads to a precise determination of sin 2 8^ at low 
momentum transfer. The running of the weak mixing an- 
gle is observed with over 6<r significance. The consistency 
of the result with the Standard Model prediction provides 
significant new limits on TeV scale physics, comparable 
in sensitivity and complementary to the best current lim- 
its from high energy colliders. 
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